Accelerating the switching of logical units by anisotropy driven magnetization 

dynamics 
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In this work the magnetization dynamics of clusters supported on non-magnetic substrates is 
shown to exhibit an unprecedented complex response when subjected to external magnetic fields. 
The field-driven magnetization reversal of small Co clusters deposited on a Cu(lll) surface has 
been studied by means of first-principles calculations and atomistic spin dynamics simulations. For 
applied fields ranging from 1 Tesla to 10 Tesla, we observe a coherent magnetization reversal with 
switching times in the range of several tenths of picoseconds to several nanoseconds, depending on 
the field strength. We find a non-monotonous dependence of the switching times with respect to the 
strength of the applied field, which we prove to have its origin in the complex magnetic anisotropy 
landscape of these low dimensional systems. This effect is shown to be stable for temperatures around 
10 K, and is possible to realize over a range of exchange interactions and anisotropy landscapes. 
Possible experimental routes to achieve this unique switching behaviour are discussed. 
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I. INTRODUCTION 

Magnetic nanoclusters deposited on surfaces present 
many intriguing properties, and have been proposed as 
building blocks for future data storage applications. The 
ability of altering the magnetization of these clusters in 
a fast and controlled way poses a great challenge. Mag- 
netization dynamics has recently been the focus of sev- 
eral theoretical and experimental investigations, both for 
bulk-like systems^— as well as nano-sized objects^— . 
From a technological point of view it is the switching of 
magnetic logical units that is of interest, since storing in- 
formation in a magnetic medium as fast and reliably as 
possible stands out as being crucial. 

Discrepancies between the properties of nano-sized ma- 
terials compared to their bulk counterpart are known for 
some time. There are several reasons why nano-sized 
objects can behave so differently compared to the bulk 
materials, e.g. due to quantum confinement effects and 
the fact that the surface to volume ratio is very large. A 
good example of this are metallic nano-particles which 
can have a completely different optical response com- 
pared to the corresponding bulk or thin-film systems. 
This is also illustrated by the famous Lycurgus cup from 
the Roman era, where gold nano-particles included in 
glass create a unique luster—. 

In many of the investigations published so far, funda- 
mental new knowledge of magnetization dynamics and 
the magnetism of nano-sized objects have been discov- 
ered. An example of this is the possibility of obtaining 
magnetization reversal on femtosecond time-scales, as re- 
ported in Ref[l|, the breakdown of the macrospin modelii 
as well as the possibility to achieve all-spin-bascd logic 
operations on an atomic levels. 

Typically, one observes one of the two known types of 
magnetization reversal: mono-domain switching or do- 
main wall motion^. In both cases, the time it takes for 
the magnetization to reverse its direction, i.e. the switch- 



ing time, is reduced when the strength of the applied field 
is increased. This is natural since a stronger field pro- 
vides a stronger driving force to reverse the magnetiza- 
tion direction. Contrary to the conventional expectation 
that a stronger applied field yields a faster switching of 
the magnetization, it is shown here that for carefully se- 
lected clusters the switching can actually be accelerated 
by decreasing the applied field. 

We present here magnetization dynamics in nano-sized 
clusters, which under special conditions deviate from 
well established connections between force, acceleration 
and speed, thus seemingly disobeying classical laws of 
physics. Our findings have no previous counterpart in 
the field of magnetism, but analogies can be drawn to so 
called non-Newtonian fluids^ (e.g. of colloids in suspen- 
sion) which demonstrate a highly non-linear response to 
an external stimulus. 



II. COMPUTATIONAL DETAILS 

In order to provide a realistic description of the sys- 
tem we have used a multi-scale and multi-code approach. 
We have first determined the ground state properties of 
these islands by means of ab initio calculations. The first 
principle study of the nanoislands' electronic and mag- 
netic structure was performed using a Green's function 
formalism within relativistic density functional theory 
(DFT), as implemented in SKKR (fully relativistic spin- 
polarized screened Korringa-Kohn-Rostocker) 2 ^ and RS- 
LMTO-ASA (real-space linear muffin-tin orbitals within 
atomic sphere approximation) 2 ^. All relativistic effects 
have been accounted for by solving the Kohn-Sham-Dirac 
equation. After having obtained ab initio site-resolved 
quantities (i.e. magnetic moments, anisotropics etc). In 
the second step, we investigate the magnetization dy- 
namics in the nanostructurcs in terms of atomisitic spin 
dynamics by means of the UppASD (Uppslaa atomistic 
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spin dynamics) package^. 

III. RESULTS 

In the present study we are particularly interested in 
the switching behaviour of the magnetization direction 
of nano-particles supported on non-magnetic substrates, 
under the influence of an external magnetic field. The 
magnetic islands are composed of 16 up to 121 Co atoms 
deposited on a Cu substrate. We have investigated the 
switching behavior of these islands in the presence of a 
static external magnetic field. The Cu substrate is a very 
suitable choice for the study of magnetic nanostructures 
supported on its surface, since it is rather inert to polar- 
ization effects due to its completely filled d-shell and its 
weak spin-orbit coupling. Hence, in a spin-Hamiltonian 

H = -- ^2 Jij rfii- rhj + ^ K l (rhi ■ e K ) 2 - B ext m t 

i¥=j ^ , J 

exchange anisotropy external 

only the Co nano-particle has to be considered, albeit 
with appropriate parameters that come from a first prin- 
ciples theory which includes also effects of the substrate. 
The systems we study consist of two atomic-layers high 
Co islands (in fec-stacking) deposited on a Cu(lll) sub- 
strate. The nanostructures are triangularly-shaped. The 
fact that the system is finite and has a low symmetry is 
reflected in the electronic structure and magnetic proper- 
ties, which has already been proven in a previous work^l. 
As a consequence, the spin and orbital moments as well 
as the exchange interactions and the magnetocrystallinc 
anisotropy have a non-uniform spatial distribution within 
the nanostructures. Since it is well known that the mag- 
netic properties depend strongly on the individual atoms' 
local environment^, we take into account the edge ef- 
fects, which have an increasing importance the smaller 
the system is. 

As a first step, we performed a thorough ab initio inves- 
tigation of the islands' electronic and magnetic structure. 
After having obtained the ab initio site-resolved quanti- 
ties, such as magnetic moments, interatomic exchange 
parameters and magnetic anisotropics, we continue with 
the investigation of the spin dynamics in these nanostruc- 
tures by means of atomistic spin dynamics simulation. 
As discussed above we do not consider the Cu substrate 
in the spin dynamics simulations, since the spin and or- 
bital polarization is negligibly small (of the order of 0.005 
fis/ atom) . Hence the tiny induced moments in the Cu 
substrate and the very weak exchange coupling with the 
Co atoms do not influence in any way the dynamics of 
the Co system. 

Using all site-resolved quantities determined from first- 
principles as initial parameters, we investigate the mag- 
netization dynamics of the system under the influence of 
external magnetic fields. The time evolution of the mag- 
netization, as described by the Landau-Lifshitz-Gilbcrt 



(LLC) equation, contains a precession and a damping 
term: 

drrii fi e ff a , _ . ^teff^ 

— — = -7 rrii x B- JJ - 7— x m, x B t J J , 
at m 

where B^ = Bi + b{ l (t) is the effective field given by 
-» OH 

Bi = — and including temperature effects intro- 
ortii 

duced by the stochastic fields b{ l (t). Since the damping- 
parameter a entering the LLC equation was not obtained 
from first-principles, its value has been varied within rea- 
sonable limits without noticeable effects on the behavior 
of the switching process. The presented results were ob- 
tained for a damping a=0.1. 

As will be shown later on, the most important quan- 
tities that determine the switching behavior are in fact 
the atom-projected magnetocrystalline anisotropy ener- 
gies (MAE). From the ab initio calculations results we 
find that the direction of the easy axes and the anisotropy 
strengths depend on the atoms positions in the cluster. 
It is possible to map the ab initio calculated anisotropy 
energies to an effective Hamiltonian where each atom 
has a unique uniaxial anisotropy. Combined, the local 
uniaxial anisotropics yield a complex anisotropy energy 
landscape of the cluster, which has a net out-of-plane 
easy axis. A depiction for the directions of the site- 
projected easy magnetization axes for the low-symmetry 
case can be seen in Fig. [JJ In addition to the direct 

(a) high-symmetry anisotropy landscape 




(b) low-symmetry anisotropy landscape 




top layer bottom layer 



FIG. 1: (Color online) Atom-projected uniaxial 
anisotropy axes: Illustrative representation for the 
different orientations of the easy magnetization axes in 
the (a) high-symmetry and (b) low-symmetry scenarios, 
for both top and bottom layers of the 16 atoms Co 
nanoisland. 

fit from the ab-inito results (Fig. [Taj, where each easy 
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axis direction points toward the corresponding global 
minima of the anisotropy energy for each specific atom, 
we have also considered a slightly altered anisotropy en- 
ergy configuration (Fig. Ilbp . This second scenario cor- 
responds to the case where for the atoms situated on 
one of the cluster's edges we chose the easy-axes to point 
in the direction of calculated local energy minima (in- 
stead of global energy minima) , in this way introducing a 
break in the anisotropy energy landscape symmetry. All 
the other atoms have their easy-axes directions pointing 
along global minima, as in the high-symmetry case. 

In this case the symmetry is lower than what was ob- 
tained from the ab-inito calculations. This latter con- 
figuration, referred in the following as the low-symmetry 
case, should be seen as an effort to make a realistic de- 
scription of symmetry-breaking effects that are likely to 
occur in these supported nano-islands, i.e. geometrical 
distortions or chemical intermixing. Before performing 
the magnetization dynamics simulations, the equilibrium 
magnetic structure is obtained by allowing the magnetic 
moments to relax in this anisotropy landscape. This re- 
sults in an essentially collinear magnetic ordering point- 
ing out of plane. For the low-symmetry case a small 
deviation from the surface normal, up to 6° in terms of 
the polar angle 9 was found. 

We investigate next the reversal mechanism of the 
magnetization, under the influence of an applied mag- 
netic field, in these magnetic nanoislands. In the fol- 
lowing we will focus our attention only on small islands 
(16 up to 121 Co atoms) since the behavior they ex- 
hibit is analogous to all sufficiently small island sizes, i.e. 
where edge effects are large and the size of the cluster 
does not allow domain formation. We study the switch- 
ing dynamics of the magnetic system under the driving 
force of external magnetic fields B ox t of different inten- 
sities, pointing along the surface normal and having an 
opposite direction to the magnetic moments' orientation, 
B xt = ( 0, 0, — B z ). During the application of the 
magnetic field, we follow the time evolution of the aver- 
age magnetization. We let the system evolve for 180 ps 
and we probe the changes in the magnetization's orien- 
tation each 100 attoseconds. Following the change in the 
orientation of the average magnetization's z-component 
we determine the switching times (t sw ) corresponding 
to different field intensities (the switching is achieved 
when the z-component of the average magnetization is 
flipped by 180° with respect to its initial orientation). 
Despite the fact that the values of the magnetic moments, 
anisotropics and exchange interactions are not uniformly 
distributed over the island, we find that a coherent mag- 
netization reversal takes place for all systems studied. 
The nanoislands behave and switch essentially as a mono- 
domain as long as the temperature is low enough. Even 
though the magnitude of the spin moments differs within 
the island, since they remain parallel during the whole 
dynamical process, one may regard them as a collection 
and, for simplicity, we refer to all these collinear spins as 
macro-spin in the following. 



First we perform the switching simulations for the 
high-symmetry case (scenario corresponding to Fig. Ila[) 
and find that as the applied field increases, the observed 
switching times decrease in a monotonic fashion for clus- 
ters exhibiting these easy-axes orientations. An example 
for the 111 Co atoms case is shown as filled squares in 
Fig. [21 This is an expected result^, as the driving force 
for the switching is stronger when the field strength in- 
creases. 

Next we consider the clusters where the symmetry of 
the magnetic anisotropy landscape is reduced (Fig. Ilbl) . 
For the 16 atoms island, the weakest magnetic field for 
which the switching occurs is ~ 1.1 T and it takes roughly 
120 ps after its application, for the island's magnetiza- 
tion to be fully reversed (open circles in Fig. This 
is considerably faster compared to the switching time 
of clusters with a symmetric magnetic anisotropy land- 
scape (filled black squares in Fig. [5]). When increasing 
the strength of the external field, the torque driving the 
switching process increases, fact that is expected to lead 
to shorter switching times. In a certain range, we find 
however that the switching times are actually shorter the 
weaker the magnetic field is. For example in the case of 
the 16 Co atoms island (empty circles in Fig [5]), for a 
field of 1.1 T the switching time is roughly 30 ps shorter 
than for B oxt =2 T. 

When increasing the applied field further (above 2 T) , 
the switching times get shorter the stronger the field is. 
For the 111 Co atoms case (empty squares in Fig. [2]), 
a 2 T applied magnetic field fully reverses the island's 
magnetization 60 ps faster than a 3 T field. For the 
larger islands, 121 Co atoms (empty diamonds in Fig. [2]), 
the switching time in an external magnetic field of 2 T is 
very short and comparable in duration to the switching 
time obtained in a much stronger field, of an intensity of 
5 T. We probed the magnetization dynamics up to very 
strong magnetic fields (up to B cxt =10 T, Fig. [5]). 

The non-monotonous dependence of the switching 
times on the external magnetic fields is quite pronounced 
for B cxt 2 T for 16 Co atoms and for B cxt ^ 3 T for 
111 Co atoms, while it becomes less pronounced for is- 
lands of 121 Co atoms (see Fig. [5]). Since for the 121 
atoms islands we map the easy-axes directions from the 
16 Co case only onto the atoms on the edges, while all 
the other atoms have the same orientation of their easy- 
magnctization axes, it is expected that the ratio of edge 
to volume effects, introduced by the anisotropybecomes 
smaller and the non-monotonous behaviour is less pro- 
nounced. Comparing the switching times of the islands 
with high-symmetric energy landscape with those of low- 
symmetric landscapes (Fig. [5]) one notices that the latter 
reverse their magnetization direction much faster. For 
relatively weak external fields, the difference is very large. 

The fact that we observe an increase in switching times 
with increasing field strength, only when the symmetry 
of the magnetic anisotropy energy is low, gives strong ev- 
idence that it is in fact the complex anisotropy landscape 
that causes this non- monotonous behavior (see Fig. [2]) . 
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external magnetic field (T) 

FIG. 2: (Color online) Switching time dependence on 

the strength of the external magnetic fields for the 
low-symmetry anisotropy landscape of Fig. Ilbl (empty 
symbols), illustrating a non- monotonous relationship 

between switching time and field strength. Filled 
squares represent the monotonous dependence of the 
switching times w.r.t field strengths, for an anisotropy 
landscape that reflects a higher symmetry of the island 

(Fig.ua. 




* (cleg) 

FIG. 3: (Color online) Trajectories of the 16 Co atoms 
macro-spin, in the magnetocrystallinc anisotropy energy 

landscape, under the influence of a 1 T 
(yellow/light-coloured line), 1.5 T (black line) and 2 T 
(red/gray-coloured line) external field, respectively. The 

color scale indicates the strength of the anisotropy 
energy (mRyd), showing that 0=90° represents the hard 
magnetization plane (i.e. the island's plane). A detailed 
illustration of the switching dynamics can be found as 
an animation in the Supplemental Material^. 



This increased resistance of the system under a stronger 
applied force represents a process analogous to the non- 
Newtonian dynamics of colloids in suspension. 

In order to elucidate this phenomenon, we now pro- 
ceed with a detailed analysis of the magnetization rever- 
sal process for clusters with a low-symmetry magnetic 
anisotropy. Since there is a coherent magnetization re- 
versal present, all the spins in the nanoisland remain par- 
allel (within very small deviations) during the switching 
process. In Fig. |31 the spin trajectories for all the atoms 
within an island would overlap, so for clarity we chose 
to represent only one spin-trajectory for each case. We 
start by investigating the change in the trajectory of the 
macro-spin, projected onto the anisotropy energy land- 
scape. We plot (in a map view) in Fig. [3l the three 
dimensional (3D) magnetic anisotropy landscape in po- 
lar coordinates, together with the paths taken by the 
magnetization vector for different field strengths. Note 
that Fig. [3] shows the 3D energy landscape in a top- view 
projection, as a function of polar angles, 9 and <fi, and 
that higher magnetic anisotropy regions are shown in yel- 
low/bright color and lower anisotropy energy regions are 
shown in purple/dark color. The changes in the polar an- 
gle 9 represent a variation in the out-of-plane component 
of the magnetization, while changes in <fi show variations 
in the in-plane component. 

The energy landscape shows a maximum in the 



anisotropy energy at coordinates 0=90° and 0=0° and 
180°, which represents the hard-magnetization region 
and corresponds to the surface plane. The fine con- 
tour lines mark equi-energy lines in the anisotropy energy 
landscape. Fig. [3] clearly shows that the energy land- 
scape is not independent of the azimuthal angle <f> which 
would be the case for a single uniaxial anisotropy of the 
nano-island. The largest energy barrier to overcome is in 
the surface plane, which is the hard-magnetization plane. 
Fig. [3] shows an important result, namely the marked 
difference in the paths that the nanoisland's macro-spin 
takes, under the influence of different strengths of the 
applied external magnetic field (illustrated for three field 
strengths: 1.0 T ycllow/light-colourcd line, 1.5 T black 
line, 2 T rcd/gray-coloured line). 

The question why the macro-spin approaches the mag- 
netically hard plane of the energy landscape quicker for 
the 1.5 T case compared to the 2 T field, can be explained 
from the trajectories shown in Fig. [3] It can be seen that 
the reason for the faster switching in the 1.5 T field is that 
this trajectory " skips" a precession and approaches much 
faster the hard magnetization plane. Translated in terms 
of the 3D anisotropy energy profile, this means that the 
resulting effective field acting on the macro-spin moves it 
to a region of the MAE landscape where the energy has 
its maximum value. Under the stronger field (e.g. 2 T), 
on the other hand, the spins follow the expected pre- 
cessional movement around the resulting effective field's 
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axis. 

For comparison, we start by describing the well-known 
case of coherent switching of a macro-spin in a single 
uniaxial anisotropy^ environment, under the influence 
of an antiparallel external field. In that case, the result- 
ing effective field acting on the macro-spin, will always 
have a constant direction along the easy magnetization 
axis. Thus, the switching will be determined solely by the 
damping torque from the effective field and the switching 
time will decrease with increasing field strengths. If the 
external field is weak, the torques generated by it will be 
counterbalanced by torques induced by the anisotropy 
field and the magnetization reversal will not occur. In 
this case, the equilibrium direction of the macro-spin will 
still be along the easy axis. 

In our system with a low symmetry of the MAE land- 
scape (Fig. llb|) . the picture is more complicated. The ap- 
plied field will still give rise to precessional and damping 
torques, but the contribution from the anisotropy field is 
significantly more complex. Due to the competition be- 
tween the different site-projected magnetic anisotropics 
in the cluster, the anisotropy fields will cause processional 
and damping torques albeit not in the same directions as 
the torques derived exclusively from the applied field. 
This can be seen very clearly in the trajectory for the 
B ext =l .0 T (yellow/light-coloured line in Fig. [3]) where 
the macro-spin makes only a small curled movement to a 
new static equilibrium position, without a magnetization 
reversal. For an intermediate field, such as B oxt =1.5 T, 
there is a delicate balance between the different torques. 
The precessional torque from the anisotropy field be- 
comes parallel to the damping torque of the applied field. 
This results in the sharp turn of the trajectory for this 
macro-spin towards the magnetically hard region of the 
anisotropy landscape. If the applied field is increased 
even more (see B ext =2.0 T in Fig. [3]), the precessional 
torque from the applied field dominates over all other 
torques. The macro-spin will, in this case, essentially be 
driven by the torque originating from the applied field, 
which dictates a precession movement. Instead of be- 
ing rapidly forced by the anisotropy torque towards the 
magnetically hard region, the macro-spin makes an ad- 
ditional revolution around the z-axis. 

Even inside the high anisotropy energy region (i.e. 
for values of 9 close to 90°), the ratio between the ap- 
plied field and the anisotropy field differs between the 
Bext=l-5 T and 2.0 T cases, even though the difference 
is not as drastic as in the early stage of the switching 
process. Both trajectories follow here roughly parallel 
paths (see Fig. [3]) but we find that the trajectory under 
the stronger field is delayed further in this region due 
to the fact that it actually crosses over the highest peak 
of the anisotropy energy landscape. On the other hand, 
the torque exerted by the lower field on the moments is 
not strong enough to overcome the highest anisotropy 
barrier. Once the moments have passed the magneti- 
cally hard region, they proceed with their precessional 
and damped motion towards the z-direction without sig- 



nificant differences between the two paths. 

A. Different anisotropy landscapes 

As we have already shown above, the shorter switching 
time with a weaker field is found for islands of different 
sizes (the largest island we took into consideration for 
our study contains 121 Co atoms). Moreover, it holds 
for different choices of the anisotropy energy landscape. 
Hence we investigated also several other cases besides 
the one depicted in Fig. [lb] and obtained similar results 
(see low-symmetry scenarios 1, 2 and 3 in Fig. [4]). All 
simulations show that a faster switching is obtained by 
decreasing the magnetic field's intensity, within a certain 
range. This is translated by the occurrence of a longer 
switching time for stronger applied magnetic fields, up 
to a certain field strength after which the switching time 
decreases with increasing fields. The only pre-requisite 
for this to happen, is the presence of a sufficiently low- 
symmetry distribution of the magnetic anisotropy over 
the island. 
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external magnetic field (T) 



FIG. 4: (Color online) Effect of weaker exchange 
interaction (by a factor 100) between the 16 Co atoms 
within the island (black stars). Different anisotropy 
energy scenarios: empty circles corresponding to the 
low-symmetry scenario of Fig. lib! empty triangles 
corresponding to random direction of the 
easy-magnetization axes. 



B. Influence of the exchange parameters' strength 

So far, we have not discussed the effect of the ex- 
change interactions on the complex magnetization dy- 
namics shown in Fig. [5] We note however, that the 
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exchange interactions obtained from our first principles 
calculations are strong enough to maintain a collinear 
arrangement, forming a macro-spin, during the stud- 
ied switching scenarios and thus their relative strengths 
are not as important for the switching as the individ- 
ual anisotropy energies. The non-Newtonian dynamic 
process is present even when the exchange parameters 
within the Co island arc different from the ab initio cal- 
culated values. In Fig. [4] we prove that the driving force 
behind the non-monotonous and accelerated switching 
behaviour is the anisotropy landscape, since this non- 
Newtonian response is recovered even in the fictitious 
case of a 100 times weaker exchange interaction between 
the Co atoms. The main difference between this case 
(black stars in Fig. |4]) and all the other considered sce- 
narios, is that here the magnetization reversal occurs in 
a non-coherent way, in the sense that the spins do not re- 
main collinear during the reversal process, but each one 
of them follows its own path under the influence of the 
external magnetic field. Thus the complex connection 
between the switching time and the strength of the ex- 
ternal magnetic field is present even when reducing the 
exchange interaction strength by a factor of 100. 




external magnetic field (T) 



FIG. 5: (Color online) Temperature dependence of the 
non-Newtonian magnetization switching, in the 121 Co 
atoms island. 



C. Temperature effects 

In Fig. O we emphasize that the accelerated magneti- 
zation switching under weak fields is present also in tem- 
perature ranges (up to 10 K) which are readily achiev- 
able in experiments. For higher temperatures, the mag- 
netization switching has a more stochastic behaviour, 
even though within certain ranges of the magnetic field's 
strength the accelerated reversal can still be obtained for 
weaker applied fields. 

Finally, we note that at sufficiently high temperature 
both the macro-spin picture and the reported switch- 
ing behavior breaks down into a more stochastic behav- 
ior, but the effect reported here is stable even above 
10 K. Both the islands size and the temperature range, 
would allow for the systems described above to be inves- 
tigated in an experimental set-up. On general grounds 
larger magnetic units, with a similar energy landscape 
as the cluster in Fig. lib! are expected to exhibit a non- 
Newtonian magnetization dynamics up to even higher 
temperatures. 



IV. CONCLUSION 

The dynamics of the magnetic clusters with a low 
symmetric energy landscape have, for a certain range 
of parameters, a dynamical response that suggests that 
well established relations between speed, acceleration and 
force seem not to apply. Stronger driving forces will in 
certain cases only slow down the dynamics. The pos- 
sibility to use this fact in technological applications is 
obvious, since faster switching of magnetic units can be 



obtained with weaker driving forces, whether it is a mag- 
netic field or a torque provided by a spin-transfer torque. 

We have illustrated a highly complex (and fast) switch- 
ing behavior for several selected cases of supported nano- 
clusters and we argue that this behavior should be 
present in a wide range of systems, even large systems 
containing thousands of atoms. Since the driving force is 
the individual atom's contribution to the total magnetic 
anisotropy landscape, which is known to be very sensi- 
tive to the local environment^, this effect is likely to be 
enhanced even further by tailoring specific clusters with 
respect to both geometry and chemical alloying. This 
can be achieved, for example, by considering clusters 
of mixed chemical composition, where local anisotropy 
axes are expected to cause an even more asymmetric 
total energy landscape of the magnetic anisotropy. An 
alternative would be to grow clusters on random alloy 
substrates, e.g. Cu x Ag(i_ x y It is known that ligand 
states from nearest neighboring atoms influence the local 
anisotropy^, hence a Co atom neighboring a Ag atom 
will have different easy axis direction than a Co neigh- 
boring a Cu atom. 

Using atomistic simulations based on density func- 
tional theory we show that this kind of non-Newtonian 
magnetization dynamics can occur provided that the en- 
ergy landscape of the magnetic anisotropy has a suffi- 
ciently low symmetry. The unexpected dynamics is ob- 
served for experimentally achievable cluster sizes (more 
than 100 atoms) and temperature ranges (above 10 K). 
We propose this effect as an enabler for faster informa- 
tion processing in technological applications, since much 
faster switching of magnetic units can be obtained with 
a weaker force, whether it is in the form of an applied 
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magnetic field or provided by a spin-transfer torque. 
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